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Abstract. We investigate an extension to the phase shift formalism for 
calculating one-loop determinants. This extension is motivated by requirements of 
the computation of Z-string quantum energies in D=3+l dimensions. A subtlety 
that seems to imply that the vacuum polarization diagram in this formalism is 
(erroneously) finite is thoroughly investigated. 



1. Introduction and motivation 

Z-strings were first discovered as solutions to the classical field equations of the 
electroweak Standard Model by Nambu [1] in the context of bound pairs of magnetic 
monopoles. Later on they were rediscovered — as independent objects in their own 
right — by Vachaspati [2] . 

The main point of interest in our study of Z-strings is their stability. If they are 
stable, they might be relevant for a variety of reasons: They would be the only static 
solitons in the Standard Model, networks of Z-strings might play an important role in 
baryogenesis [3], and they would be a possible source of the primordial magnetic field. 
For a general overview of applications and properties of Z-strings along with a large 
collection of references cf. [4]. 

Because there are no configuration with nontrivial topology in the electroweak 
model, stability of Z-strings is only possible on energetic grounds. Classically, Z-strings 
are unstable for the physical values of the parameters of the electroweak model. Since 
fermions bind strongly to the core of the Z-string, it may be possible to achieve stability 
by occupying N fermionic bound states that are generated in the background of the 
Z-string so that the the resulting total energy of the system is less than the energy of 
TV free fermions. A consistent fi expansion then requires to include the contribution 
of the fermionic determinant to the energy. 

In D=3+l dimensions, renormalization issues make investigation of the Z- 
string difficult since it is non-perturbative. Investigations have either failed to draw 
convincing conclusions [5] or focused on D=2+l dimensions [6l [7] . Here we will discuss 
attempts to solve the problems posed by these earlier investigations. 
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The structure of this paper is as follows: In section [2J we discuss our method 
for computing the fermion determinant and the problems posed by the Z-string for 
our approach. We discuss the possible solutions in D=2+l dimensions and the 
extension necessary for D=3+l dimensions. In section [3] we investigate whether 
the extension proposed for the Z-string works in the case of a simple boson model. 
Problems involving the second order Born approximations are discussed and resolved 
by comparing different formulations of our approach. In section [4] we present some 
conclusions and an outlook. 

2. Technical prerequisites 

The Z-string has the following structure: 

*- (* )-( (1 » 

with </> the iso-spinor Higgs field and Z the only non-vanishing component of the 
SU (2) x U{1) gauge fields. As usual, p and <p denote the radial and angular coordinate 
in the plane perpendicular to the Z-string and v is the vacuum expectation value 
of the Higgs field. The string configuration is defined in two spatial dimensions. 
It is translationally invariant in any additional space coordinate, which we call flat 
dimension(s). Finite classical energy (for D=2+l, or classical energy per unit length 
in D=3+l) requires both /# — * and fa — > for p — > 0. Furthermore, the kinetic 
energy of the Higgs must fall off sufficiently fast at large distances, 

Kfy - i 9 Z^\ 2 = 0(l/p 2+e ), e > for p oo. (2) 

In particular, this implies that <j^<j) = \ Jh\ 2 — > 1 and fa — > 1 at p — > oo. The condition 
([2]) clearly mixes different orders in perturbation theory (PT) and already indicates 
that if one sticks to a fixed order in PT one will end up with IR divergent quantities. 
In the systematic expansion of arbitrarily many fermion species and at next to leading 
order in h, the quantum correction to the energy, i.e. the vacuum polarization energy, 
may be computed from the fermion determinant in the background of the potential 
V(x) that is generated by the string. In turn, this determinant is given in terms of 
scattering data from this potential [8j . This formulation is effective for renormalization 
because the n*' 1 order contribution in the Born expansion equals the Feynman diagram 
for the fermion loop with n insertions of V. To render the integral over scattering data 
finite, we subtract enough Born terms and add them back in as Feynman diagrams. 
The renormalization of the latter is adopted from perturbation theory and is standard. 
Hence the vacuum polarization energy is given by 

1 1 r N 

= ^ E fM'-) - ^ / dfc E [m*o - E *M (*) 

b.s. J M n=l 

N 

+ Y i E% + E C i, (3) 

71=1 

where w^ s - denotes the fermion bound state energies, 8m the (full) phase shift in 

angular momentum channel M, 5^ its rt th Born approximation, N denotes the 
number of Born approximations necessary to render the momentum integral finite, 

(n) 

Epf, is the energy contribution computed from Feynman diagrams with n external 
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legs and i?cT is the energy resulting from the counter terms. The functions fx 2 vary 
with the number of flat dimensions [9 . For the Z-string in D=2+l, the number 
of flat dimensions is and we obtain = m — \u>\, f2(k) = — fc/v k 2 + to 2 . 

For D=3+l, the number of flat dimensions is 1 and we obtain fi(u) = (w 2 (l + 
ln(o; 2 /TO 2 )) — to 2 )/(47t), / 2 (fc) = (k/2n) In (1 + k 2 /m 2 ). Of course, the Feynman 
diagram and counter term contribution also differ for D=2+l and D=3+l, but these 
expressions are well-known and will not be repeated here. When computing the Born 
approximation for scattering in the Z-string background, one encounters two problems: 
The full scattering problem is decomposed into distinct channels of finite size. The 
decomposition of the full scattering problem utilizes the generalized axial symmetry 
of the Hamiltonian h in a Z string background: 

[h, J 3 - nT 3 7 5 ] - (4) 

where J3 is the fermionic angular momentum, T3 is an isospin generator and n is the 
(integer) flux of the Z-string. Since 75 does not commute with the free Hamiltonian ho, 
we have [ho, J3 — nT 3 j 5 ] 7^ 0. Thus, in contrast to standard problems, the symmetries 
of the problem with background field are not a subset of the symmetries of the free 
problem, but instead are incompatible with them. Thus, the setup for perturbation 
theory (which considers the Hamiltonian h f = ho + eV) has no axial symmetry at all. 
The second problem is related to the aforementioned IR problems of strings. Only 
special (gauge invariant) combinations of Feynman diagrams of different perturbative 
orders result in IR finite quantities. However, Born approximations correspond to 
sums of Feynman diagrams of definite order in PT and hence most likely will suffer 
from IR divergences. 

The case of D=2+l is special. The only counterterm that has a divergent co- 
efficient — 1 4>\ 2 — is of a definite order in PT. In this case we are able to use gauge 
invariance to argue that any background that has the same \<f>\ 2 will suffer from the 
same divergence and hence the Born approximation for some other fake background 
can be used to subtract the large-momentum tails of the full phase shifts. So we can 
choose a background that has J3 as a symmetry generator but no IR divergences [7]. 

In D=3+l this procedure does not work because we now have more counter terms 
with divergent coefficients that mix different orders of PT. Furthermore, since there are 
additional constraints on the profiles, it is not clear how to construct fake background 
solutions which have no winding and thus the same symmetry as the free Hamiltonian 
- if it's possible at all. 

At this point it is worthwhile to remember why we want to subtract Born 
approximations in the first place: First, we have to subtract off the large-momentum 
behaviour of the summed phase shifts in order to ensure UV convergence of the 
momentum integral. Second, with Born approximations we know exactly what to 
add back in so that the overall value of the determinant does not change, namely the 
corresponding Feynman diagrams. Keeping these reasons in mind, we can extend 
our formalism substantially: Instead of subtracting the Born approximations of 
the original theory with the original background fields, we can subtract off Born 
approximations from an arbitrary theory as long as a) we add back in the Feynman 
diagrams from the same theory and b) that theory has the same divergences (e.g., in 
dimensional regularization) as the original theory. In the remainder of this paper, we 
test this proposition. 
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3. Fake subtractions 

We consider two bosonic theories in D=3+l that differ in the string- like background 
potential they're coupled to: 

1 m 2 1 

A,2 = -d^ d»cf> ~ — <P 2 ~ -C71,2 <f (5) 

with o~i(x) = a.\ exp (— p 2 /w 2 ), 02(2?) = at2 exp (— pjw^). Demanding identical UV 
divergences requires 

J d 3 xcri(x) = J d 3 xa 2 (x) and J d 3 xa 2 (x) = J d 3 xaj(x). (6) 

If our proposal of a fake Born subtraction works, the renormalized vacuum polarization 
energy in ^ should be the same, whether we use the original model 1 or its fake version 
2 for the subtraction. Since the bound states are not affected by the subtraction, it is 
sufficient to consider the part 

El = Udkk\n(l + ^)j2l [Sm - £ S$ (fc; <7,)]+X; 4t re „(fc; °i) (7) 

J ^ ' M n=l n=l 

where the full phase shift is always computed in the original model 1, while the 
subtraction is made with model 1 Born approximations (E\) or with the fake model 
2 Born approximations (E2), respectively. Our proposal amounts to claiming that 
Ei = Ei. 

It can be shown, using a Bessel function identity [8], that X^m=o oc 
{o~i)k q ~ 2 where q is the number of non-ffat spatial dimensions. By numerical 
calculation we can furthermore show that 

\ J dkkln^l + ^Y.l^M^^-S^ik;^)} =41, e „(^)"41, e „(^i)- 

(8) 

This shows that as long as the renormalized determinant is finite, both ways of 
subtracting Born approximations lead to the same result. In the left panel of figure [1] 
we see that the renormalized determinant is going to be finite, since there is no large- 
momentum tail that might impede the existence of the momentum integral. The 
right panel, however, shows something troubling: Ordinarily, we would think that ([7]) 
with just the first term of the Born series subtracted still contains the log-divergent 
2 nd order Feynman diagram. The Feynman series then suggests the same large k 
behaviour for the integrand of {7} when the expression in square brackets is replaced 
by the second order term of the Born series. From the right panel of figure [TJ one 
can clearly see that this expectation is wrong: both contributions to the integrand 
in {2} fall off faster than 1/k 3 individually. At first sight this indicates that the 
resulting integral would be finite [101 E] • This is in strong contradiction to the fact 
the equivalent 2™ d nd order Feynman diagram is indeed UV divergent. As we will 
demonstrate in the rest of this paper the catch comes from an incorrect treatment of 
the non-uniformly convergent sum and integral in ([7]). 

First of all, it is worth noting that the second order Born approximation phase 
shifts are not small by themselves, but they sum up to something small, cf. figure 
[2] Secondly, it is worthwhile to recall the formulation of functional determinants in 
terms of scattering data P~2J [T3] . Strictly this is possible only in the upper half of the 
complex k plane and not on the real k axis. Furthermore, that derivation also suggests 
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Figure 1. The left panel shows k In (l + fc 2 /m 2 ) i [<5 M (>; o"i) -<5^ (fc; <ti) -5^ (k; (Tj)] fori: 

1, 2. The right panel shows the decomposition furthermore into k In (l + fc 2 /m 2 ) ^ M ^ [<5m ""l) 
«&Wl)] and fe In (l + fc 2 /m 2 ) £ M £ 5$ (*i <r<)- 




Upper channel limit Upper channel limit divided by k 



Figure 2. The left panel shows S^(k) = (k; (72) for different values of k, the right panel 

shows k 2 sP(k) plotted vs. iV/fc. 



that instead of integrating over a channel sum, one should integrate over momentum 
first to get a per-channel contribution to the energy and then sum over channels. 
The analytic properties of scattering data ensure that the contribution of a prescribed 
channel to the vacuum polarization energy can be identically computed as an integral 
over real or imaginary momentum. We have verified this result numerically. Of course, 
the corresponding integrands are not expected to be identical. Consequently, the 
momentum integrands obtained from summing over channels first are expected to be 
different as well. The momentum integrand of the second order Born approximation 
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Momentum Upper channel limit 

Figure 3. The left panel shows t 5^ (t; <Tj) (which is the expression on the imaginary axis 
corresponding to k In + fc 2 /m 2 ) 5^- on the real k axis) for different values of t and i = 1,2. 

The right panel shows y"]j^_„ Jffi (t; 1T2) plotted vs. iV for different values oft. 



phase shift§| on the imaginary axis shows exactly the expected 1/t falloffjj] for a 
logarithmically divergent integral, cf. left panel of figure [3] In the right panel, one can 
see that the 2 nd order Born approximation on the imaginary axis have the same sign 
for all channels and no cancellations are present. Moreover, our general experience 
together with our numerical investigations in this particular case suggest that the 
sum over channels is uniformly convergent and we can freely interchange the order of 
channel summation and momentum integral. 

The left panel of figure 0] shows clearly that the problem on the real axis is the 
interchange of integration and summation, because by integrating first over momentum 
and then summing over channels we obtain the expected logarithmic divergence, 
whereas in the other order the result is apparently finite. The right panel shows 
how the (identical) values on the real and imaginary axis are obtained — on the 
imaginary axis we see a monotonic approach to the final value, while on the real axis 
we see a zero crossing. This zero is not accidental, but rather a consequence of the 
fundamental sum rule 

J dkk5^\k) = 0, (9) 

proved in |14j . On the imaginary axis no corresponding sum rule exists. The 
consequence of this zero crossing can be seen in figure El From the right panel of this 
figure one can see clearly that while for zero flat dimensions the dominant contribution 
to the per-channel contribution of the energy comes from the low-k region, for one 
flat dimension the dominant contribution comes from the region k ■ w > M, where w 
is the characteristic width of the potential. This result explains why the interchange 
of summation and integration fails on the real axis for one flat dimension (and not 
for zero); if one integrates first over momentum, one obviously is able to catch the 

§ More precisely one should talk about the imaginary parts of the logarithm of the Jost function, 
because the concept of a phase shift cannot be extended to the complex k plane. 
| The momentum on the imaginary axis is denoted t, i.e., k = it with t real. 
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Figure 4. The left panel shows J dk k In (1 + fc 2 /m 2 )5^ (A;) - denoted as 'real axis' - and J dt 1 6^ (t) 
- denoted as 'imag. axis' - for both <j\ and (T2 as potentials. Note that for larger values of M, the curves 
for <ri 2 coincide and show the expected l/M fall off. The right panel shows the momentum integrand 
for both real and imaginary axis for channel M = 10. 
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Figure 5. The left panel shows another scaling rule; the most important point is that the zero crossing 
moves outwards on the real k axis as M increases. The right panel compares the integrands of the 
per-channel contribution to the energy with zero and one flat dimension to the integrand of the sum rule. 



dominant large momentum contribution. If one sums over channels first, however, 
one finds that for M > k ■ w the phase shifts drop off exponentially. Hence it seems 
sufficient to sum up to M max w (2...3)fc • w. Since the momentum integral will 
be terminated at some finite value fc max , this procedure obviously misses important 
contributions from channels with M > M max . Hence, integration and summation 
can not be interchange Another piece of evidence is that the momentum integrand 

K During our investigation of the second order Born approximation phase shifts, we noted that 
k s S\y(k) is not a function of k and M individually, but is just a function of k/M (plus corrections, 
but those are very small for M > 1). This allows to map M W v ' a ^ e (l eacn ng order) Euler- 

Maclarin formula to f dk k sff (k) which might ultimately explain why on the real axis the sum over 
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obtained from summing over channels first is — as far as we can tell from our numerical 
investigations — strictly positive for any prescribed finite momentum, which is in 
direct contradiction to the sum rule from 

4. Conclusions 

We have shown that for finite quantities one can successfully replace Born 
approximations in one theory by Born approximations from another as long as the 
Feynman diagrams show the same divergences. This is a major improvement over 
existing formalisms: One does not have to deal with problems that originate from 
non-vanishing structures at spatial infinity, but still one knows exactly what to add 
in again to implement renormalization conditions known from perturbation theory. 
In these conference proceedings, we have focused on a pathological case where the 
ultimate object of interest does not strictly speaking exist (as it is UV divergent); 
but the apparent convergence of the vacuum polarization diagram in the phase shift 
formulation was important to understand. This (erroneous) convergence only occurs 
when formulating the vacuum polarization energy in terms of scattering data for real 
momenta. Analytically continuing to the imaginary axis yields the UV structure that 
is consistent with the analysis of the Feynman diagrams without any further subtleties 
like ordering of limits. Though this favors the formulation in terms of imaginary 
momenta it contains drawbacks that we did not go into — one needs to sum a lot 
more channels and it may pose challenges for fermions. 
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